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' Motivated by the recent measurements of non- leptonic B" decays by CDF and LHCb 

collaborations, especially the large B{B® — > 7r + 7r - ), we revisit the hard spectator and anni- 
hilation corrections in B® decays within QCD factorization approach with two schemes for 
the possible parameters for the known end-point divergence appeared in the estimation of 
(T) | the hard spectator and annihilation diagrams. The first one is to conservatively estimate 

the possible contributions by parameterization (scheme I); another one uses an infrared 
| finite gluon propagator (scheme II) to regulate the end-point singularity. In scheme I, with 

the constraints from the measured B® — > PP (VV) decays, two (four) restricted solutions 
of the parameters spaces are found. In scheme II, we find that most of the theoretical pre- 
dictions agree well with the experimental data with single parameter m g ~ 0.5GeV. How- 
! ever, within both schemes, B(B® — > (p<p) are always much larger than B(B® — > K*°K*°) 



in contrast to the experimental results B(B® — > <p4>) ~ B(B® — > K*°K*°). It is noted that 
the pattern B(B® — > (fxfS) > B(B® — > also persists in other theoretical framework, 

thus the present experimental result B(B® — > cpcp) ~ B(B® — > K*°K*°) rises a challenge 
to theoretical approaches for B non-leptonic decays. Further refined measurements and 
theoretical studies are called for to resolve such a possible anomaly. 
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1 Introduction 



The pure annihilation non-leptonic B meson decays, without the interference induced by spec- 
tator diagrams are very suitable for probing the strength of annihilation contribution and 
exploiting the related mechanism. Recently, CDF and LHCb Collaborations have reported the 
evidence of the pure annihilation decay — > ty + tt~, with a significance of 3.7a and 5.3cr, 
respectively, 

B(B° S 7T+7T-) = (0.57±0.15(stat) ±0.10(syst)) x 10" 6 , CDF [TJ (1) 
B(B° S tt+tT) = (0.98^ (stat) ±0.11(syst)) x 1(T 6 , LHCb [2] (2) 

as well as the branching fraction of the pure annihilation decay B d — > K + K~ , 

B(B° d ^K + K~) = (0.23±0.10(stat) ±0.10(syst)) x 10~ 6 , CDF [TJ (3) 
B(B° d ^K + K~) = (0.13l{j;gjj(stat) ± 0.07(syst)) x 10~ 6 . LHCb [2] (4) 

Averaging the experimental data Eqs. ([T|) and ([2]) , Heavy Flavor Averaging Group (HFAG) 
gives 

B(B° 8 tt+tT) = (0.73 ± 0.14) x 10" 6 . HFAG [3] (5) 
Averaging the experimental data Eqs. ([3]) and (j4j) roughly, we get 

B(B° d K + K~) = (0.16 ±0.08) x 10~ 6 . (6) 

Theoretically, the pure annihilation non-leptonic B meson decays are expected much rare 
with a branching fractions at the level 10~ 7 or less due to the fact that the annihilation 
corrections are formally AQCD/f^b power suppressed. While, together with the chirally en- 
hanced power corrections, they offer interesting probes for the dynamical mechanism governing 
these decays and exploration of CP violations, and therefore attract much more attention 
recently [U EJ El E]- Unfortunately, in collinear factorization approach, the calculation of 
annihilation corrections always suffers from end-point divergence. Within the perturbative 
QCD (pQCD) approach [8], such divergence is regulated by introducing the parton transverse 
momentum kj- at expense of modeling additional kx dependence of meson distribution func- 
tions, and present a large complex annihilation corrections [6j |8]. The most recent renewed 
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pQCD estimations of B(B° — > tt + tt^) and B(B d — > K + K~) [7] are in good agreement with 
the CDF and LHCb measurements, however, a systematic examination combined with other 
correlated decays in the same framework is not available yet. In the soft-collinear effective 
theory (SCET) [S], the annihilation diagrams are factorable and real [TU] to the leading power 
of 0(a s (m b )A QC D/m b ). 

In the QCD factorization approach (QCDF) [11], there are mainly two ways to deal with 
the end-point singularity in weak annihilation calculation: (i) scheme I, parameterization in 
a model independent way [12] with at least two phenomenological parameters introduced, for 
example Xa = Jq dy/y = ln(mb/A/ l )(l + pa^ a )\ ( n ) scheme II, using the infrared finite gluon 
propagator p31 [14], for example 1/k 2 -»■ l/(k 2 - M g (k 2 ) + ie). 

As a popular way, the scheme I is widely used in the theoretical calculations [121 EH EES! E] . 
Fitting to the data of B u ^ — > PP decays, a favored parameter value choice "Scenario S4" is 
obtained in Ref. [12]: p^ d (PP) ~ 1 and <p^ d (PP) ~ —55°, which leads to the prediction 

B(B° d K + K~) = 0.070 x 10~ 6 . (7) 

Assuming the default values of pa{PP) and 4>a{PP) in B s decays are similar to that in B u ^ 
decays, Cheng et al. give the prediction [T7] 

B(B» -+ vr+O = (0.26l°;Ko£) * 10" 6 • (8) 

It is noted that above QCDF predictions are significantly smaller than the measurements 
Eqs. ()5]) and Especially, the default value B(B° S — > 7r + 7r~) = 0.26 x 10~ 6 is about 3.4a lower 
than the experimental data (0.73 ±0.14) x 10~ 6 , which implies possible much larger annihilation 
contributions in B s decays than previous prospect. Using the CDF results in Eq. ([T]) solely, 
a detail study about such topic has been performed by Zhu [4]. Assuming universal values of 
Pa{PP) and 4>a{PP) for B^ and B s decays, it is found that QCDF is hardly to provide results 
in agreement with all of the well measured B — > PP decays. Then, if the recent measurement 
of LHCb in Eq. ([2]) is considered, the tension would be further enlarged, which may imply the 
parameters pa and <pA are non-universal in Bd and B s decays. So, it is worthy to fit their values 
with available data of B^ and B s decays, respectively, and update the QCDF predictions. 

Within the scheme II, the formula of annihilation corrections for B — > PP and PV decays 
have been given in Ref. [14]. In this scheme, with the only one input parameter effective 
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gluon mass scale m g = 0.50 ± 0.05GeV, the theoretical predictions of the observables for 
B U) d — > nK,pK and irK* decays are consistent with the experimental data [H]. So, it is 
deserved to check if its predictions for the pure annihilation decays are in agreement with the 
same effective gluon mass scale parameter. Furthermore, the pure annihilation P^,s — > VV 
decays, which involve more observables, may play an important role to test the methods of the 
end-point singularity regulation. So, in this paper, we calculate the annihilation corrections 
related to B^s VV decays with the infrared finite gluon propagator. 

In Section 2, we briefly review the annihilation contributions within QCDF. In Sections 3 
and 4, with schemes I and II for the end-point divergence regulation, we revisit B Q S — > PP , PV 
and VV decay modes, respectively. In our evaluations, the pure annihilation B s non-leptonic 
decays and the related well measured ones are examined simultaneously. Section 5 contains 
our conclusions. Some amplitudes of B® decays and the theoretical input parameters are 
summarized in Appendix A and B, respectively. 



2 Brief review of the annihilation corrections within QCDF 

In the Standard Model (SM), the effective weak Hamiltonian responsible for b — > p transitions 
is given as [T8] 

10 



v ub v: p (do? + c 2 o u 2 ) + v cb v; p (dot + c 2 o c 2 ) - v tb v* ( 



i=3 



-h.c, (9) 



where V qb V* s (q = u,c and t) are products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements, Ci the Wilson coefficients, and Oi the relevant four-quark operators . 

With the effective weak Hamiltonian Eq. (Q, the QCDF has been fully developed and ex- 
tensively employed to calculate the hadronic B meson decays. The basic theoretical framework 
of B uAs PP, PV and VV decays could be found in Refs. PH [121 [15l HQ E] . In this paper, 
we adopt the same convention and formula given in Refs. [121 HSj , except for some corrections 
pointed out by Ref. |19j . It is noted that the strength and associated strong-interaction phase of 
annihilation corrections and hard-spectator scattering contributions are numerically important 
to evaluate the branching ratios, the CP asymmetry and the polarization observables. Unfortu- 
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nately, such power correction terms always suffer from the endpoint divergences, which violate 

the factorization. To probe their possible effects conservatively, the endpoint divergent integrals 

are treated as signs of infrared sensitive contribution and usually parameterized by [121 f!5] - 

1 dx ^ f 1 , lnx I.,, , 9 f 1 dx ^ . , 

~^X A , / dx—^—(X A )\ / - 2 ^X L , (10) 

X Jq X z Jo x 

X A = (l + pe^)ln^, X L = (l + pe^ (11) 

with Aft being a typical scale of order 0.5GeV, and p, <fi being unknown real parameters. Xjj is 
treated in the same manner. The different choices of the parameters space of p and correspond 
to various scenarios, which have been thoroughly discussed in Refs. [121 HS1 HS1 ITT] . 

Fitting the fruitful experimental measurements of B u < i — y PP, PV and VP decays, a favored 
scenarios S4 is obtained in Ref. [T2]. Furthermore, the fitted p and <fi for B u d — > decays 
are also given in Ref. [T5l IT6] . Their results are summarized as 

p PP = l, tf P = -55°; (12) 
P7 = l, 07 = -20°, P r = l, 0r = "7O°; (13) 
pf *-< K * R * = 0.78 , *> M * = -43° ; p)f = 0.65 , = -53° . (14) 

Assuming the default values of p^ and 0^ in the B s decays are similar to that in B u ^ decays, 
Ref. [T7j takes the values 

pf P = l, 0f P = -55°; (15) 
pf y = 0.85, ^ = -30°, P r = 0.9, tf p = -65°; (16) 
pr = 0.70, 0r = -55°, (17) 

as the inputs for the B s decays. In this paper, we denote above parameter space as "scenarios 
54" for convenience. It is noted that some non-leptonic B° s decays have been well measured in 
recent years, such as B° s — y 7r + 7r~, tt~ K + ', K~K + , K*°K*° and (fxfi decays. So, it is worth to 
check above parameter values and refit them with the updated data of B s decays. Furthermore, 
without the interference induced by spectator diagrams, the pure annihilation non-leptonic B s 
meson decays, such as B° s — y rnr, pit and pp decays, are very suitable for probing the strength 
of the annihilation corrections and related mechanism. So, in this paper, we mainly pay our 
attention to such two types of B s decays. 



Table 1: The numerical results for the branching fractions [xlO 6 ] and the direct CP violations 
[xlO -2 ] of Bg — > txtx , pir , pp , ir~K + and K~K + decays in each scenarios. 
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3 B° s ->• PP and PI/ decay modes 
3.1 Within Scheme I 

With the annihilation parameters of scenarios S4 for B° s — > PP and PV decays given by 
Eqs. (fT5|) and (jl~6l) . and the other input parameters listed in Appendix B, the predictions 
for the observables of pure annihilation decays B° s — > tttt , prr and the well measured decays 
B° s — > 7r~K + , K~K + are given in the third column of Table [TJ The theoretical uncertainties 
are mainly induced by the three parts: quark masses , CKM elements and decay constants, 
form factors. We first scan randomly the points in the allowed ranged of the input parameters 
of the three parts, respectively, and then add errors in quadrature. 

Our QCDF results of scenarios SA listed in Tables Q] are consistent with the former predic- 
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1: The dependence of B{B° S -> vr+vr" , K + K~ , ir~K + ) and A C p{B° s -»• vr - !^) on the 
5^ PP with different /9 PP values. The dashed lines correspond to the error bars (1.68a). 



tions given in Refs. [121 E] > an d the small difference is mainly induced by the different input 
values and some corrections mentioned above. One may find most of the predictions agree 
well with the experimental measurements. However, we again find the theoretical prediction 
B(B® — > 7r + 7r _ ) ~ 0.21 x 10~ 6 , which agrees well with the other theoretical results such as 
~ 0.26 x 10" 6 in Ref. pi] and ~ 0.155 x 10" 6 in Ref. [12], is about 3.7a lower than the average 
of experimental data (0.73 ± 0.14) x 10" 6 . 

The Fig. [TJ shows the dependence of the measured observables of B® — > PP decays on the 
phase 0f p with different pf p values. From Fig. [T](a), one may easily find that the annihilation 
correction with the nominal annihilation parameter value pf p ~ 1 is hardly to account for 
the measured large B(B° S — > tt + it^) within errors, and a larger pf p is required. For the other 
measured observables, as Figs. [T] (b), (c) and (d) show, the large annihilation correction is 
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not essential, and is allowed. So, it is worthy to evaluate the exact values of the annihilation 
parameters with the constraints from the available experimental information of B® — > PP 
decays. 
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Figure 2: The allowed regions for the annihilation parameters cj) pp and p pp under the con- 
straints from the observables labeled in figures, respectively, (Figs, (a) and (b)) and their com- 
bination (Fig. (c)). 



To keep the predictive power of the QCDF framework, we assume the annihilation param- 
eters are universal for all of the B° s — > PP decay channels in this paper. Under the constraints 
from B(B° S — > 7T + n~ ,tt~ K + , K + K~), Acp{B® — > ir~K + ) and their combination, the allowed 
regions for the annihilation parameters (j) PP and p PP are shown by Figs. [2] (a) and (c), re- 
spectively. From Fig. [2] (a), we find that the traditional treatment p pp ^ 1 is allowed by 
the experimental results of B(B° S — > tt~K + ,K + K~) and Acp(B® — > ir~K + ), but obviously 
excluded by recent experimental measurements B(B Q S — > tt + 7t~) = 0.73 ± 0.14. Combining 
the constraints from above four observables, as Fig. (c) shows, the annihilation parameters are 
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tightly restricted, and two solutions, named S PP A and B, respectively, are obtained as Q 

pf p = 2.5 ± 0.8 , <f a p = -84° ± 21° ; (S PP A) 
p PP = 3.5 ± 0.3 , <f>? p = 1 16° ± 9° . (S PP B) 
Both of them imply a large annihilation correction is essential to accommodate the measured 
B s — > PP decays. 

As a comparison, we also evaluate the values of the annihilation parameters in B® — > 
PP decays with the constraints from the well measured B^ — > ir + K~ and recent measured 
B® — > K + K~ decays. From Fig. [2] (b), we find <p pp and p pp are bounded strongly by the 
precise experimental data of the branching fraction and direct CP asymmetry of B% — > tt + K~ 
decay. While, the constraint from B(B^ — > K + K~) is weak due to the rough measurement. 
Corresponding to the allowed region for (p pp and p pp shown by Fig. [2] (c) , we get the numerical 
results 

p pp = 1.2±0.3, </> pp = -48°±16°, (19) 

which is similar to the result of scenario SA given by Eq. ( fl"2"]) . while significantly different from 
those by B s decays <p PP and p PP in Eq. (TT5T) . 

For B° s — > PV decay modes, so far, there is no available experimental measurement could 
be used to put a constraint on the annihilation parameters therein. Thus, in our numerical 
evaluations, we assume that the differences between p(4>) PV ' VP and p(4>) pp in B° s decays are 
the same as the one in B^ decays of scenario ^4 given by Eqs. fTl~2T) and ([TBI , i.e., 

P 7 = PV = P P S P , 0r = 0f P ±35°, 0r = 0f P -15°. (20) 

With the default values of p pp < PV < VP and (f) PP > PV V p given by Eqs. <^> and ([20]) as inputs, 
we present our results of the observables in fourth and fifth columns of Table [TJ We find that 
B(B® — > 7T + 7r _ ) could be enhanced to 0.7 x 10~ 6 to match the experimental data with large 
annihilation parameters within QCDF. Furthermore, all of the other theoretical results are in 
good agreement with the experimental data. The branching fractions of B° s — > PV decays, 
which may play an important role to confirm or refute the large annihilation corrections, are 
too small ~ O(10~ 8 ) to be measured very soon at LHCb. 



x Out fitting for the annihilation parameters is performed with that the experimental data arc allowed within 
their respectively 1.68<7 (~ 90% C.L.) error bars, while the theoretical uncertainties are also considered by 
varying the input parameters within their respective regions specified in Appendix B. 
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3.2 Within Scheme II 

Alternative to the way of the parameterization, the end-point divergency could be regulated 
by an infrared finite dynamical gluon propagator, which has been successfully applied to the B 
meson non-leptonic decays [T3"j [T4"l |2~U] . In this paper we adopt the Cornwall's description for 
the gluon propagator, which is given by (21] 

with the dynamical gluon mass 

' q 2 +im 2 a \ -. _ 12 



M 2 (q 2 )=m 2 



In 



1? 

QCD 



ln'^ 



A 2 
QCD 



(22) 



where q 2 is the gluon momentum square, which is space-like in the hard spectator scattering 
contributions and time-like in the annihilation corrections. The corresponding strong coupling 
constant reads 

4:71 



/3 ln 



"a 5 

^QCD 



where /3q = 11 — ^rif is the first coefficient of the beta function, and n/ the number of active 
flavors. 

With the description given above, the amplitudes of the hard spectator scattering contri- 
butions and annihilation corrections in the B — >■ PP and PV decays have been derived in 
Ref. [H]. Within this scheme, it is found that the hard spectator scattering contributions are 
real and the annihilation corrections are complex with a large imaginary part [H] . Moreover, 
the strength of the annihilation correction is sensitive to the sole input parameter, the effec- 
tive gluon mass scale m g , which typical value is 0.5 ± 0.2 GeV obtained by relating the gluon 
mass to the gluon condensate [21] . Interestingly, in B meson system, with the constraints from 
B u ,d , nK* and pK decays, a reasonable similar result m g = 0.5 ± 0.05 GeV is found in 

Ref. [H]. So, as a crossing check, it is worthy to evaluate the value of the effective gluon mass 
scale m g in B° s decays. 

With m g = [0.3, 0.7] GeV allowed, the dependences of the measured observables B(B° S — > 
7r + 7r~ , K + K~ , tt~K + ) and Acp{B® — > n~K + ) on m g are shown in Fig. [3J From Fig. E] (a), we 
find that a small m g ~ 0.43 GeV, which would lead to large annihilation corrections, is required 
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(c) (d) 

Figure 3: The dependence of the measured observables of the B® — > PP decays on the effective 
gluon mass scale. The dashed lines correspond to the error bars (1.68cr). 



by large experimental data B(B° S — > 7r + 7r~) = 0.73 ±0.14. While, as Fig. |3](b) shows, a relative 
large m g ~ 0.52 GeV could result in a good agreement between the theoretical prediction and 
experimental data for B(B° S — > K + K~). With the experimental errors considered, one also 
could find a common solution 

m g = 0.48 ±0.02 GeV. (24) 

Where, its upper limit is dominated by B(B° S — > n + 7T~), and the lower one is dominated 
by B(B° S — > K + K~). Moreover, due to that a larger annihilation strength is required by 
B(B® — > 7r + 7r~), such a solution is a bit smaller than the finding m g = 0.5 ± 0.05 GeV in B u d 
decays [14] . although they are also in agreement. Due to the dominance of the tree contribution 
«i in the amplitude of B® — > tt + K~ decay, the effect of the annihilation contributions related 
to m g to B{B® — > tt + K~) is negligible, which can be seen from Fig. |3] (c). Furthermore, as 
Fig. El(d) shows, because of the large experimental error, the constraint from Acp(B® — > n + K~) 
on m g is weak too. 

Taking m g = 0.48 GeV, our numerical results for the observables are listed in the sixth 
column of Table [TJ One may find all of the results are in good agreement with the experimental 
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data, and most of them are similar to the scenarios S A and S B in scheme I. Within 
scheme I, Acp(B° s — > n + p~) and Acp(Bg —> p + ) could be large due to the assumption 
that (j)f v ^ <pY F ', which can be seen from Eq. (1201) . However, within scheme II, we find that 
Acp{B° s — > 7r + p~) = Acp{B® — > 7r _ p + ) w with any value of m g . The significantly different 
predictions for such two observables within two schemes will possibly be judged by the running 
LHCb and upcoming super-B experiments. 

4 — > VV decay modes 

Compared with B — > PP and PV decays, the B — > VV decays involve more observables, which 
are sensitive probes for testing the SM and various calculation approaches. Recently, B® — > 
K*°K*° and 00 decays have been measured by LHCb and CDF collaborations, respectively, 

B(B° S -> K*°K*°) = (28.1 ± 4.6(stat.) ± 4.5(syst.) ± 3A(f s /f d )) x 1(T 6 , 
< f L (B° s K*°K*°) = 0.31 ± 0.12(stat.) ± 0.04(syst.) , LHCb [22] (25) 
f ± (B° s K*°K*°) = 0.38 ± 0.11(stat.) ± 0.04(syst.) ; 

and 

B(B° S 00) = (23.2 ± 1.8(stat.) ± 8.2(syst.)) x 10~ 6 , 
< f L (B° s 00) = 0.348 ± 0.041(stat.) ± 0.021(syst.) , CDF j23] (26) 
f ± (B° s 00) = 0.365 ± 0.044(stat.) ± 0.027(syst.) . 

Because of the left-handedness of the weak interaction and the fact that the high-energy QCD 
interactions conserve helicity, the SM expects the dominance of the longitudinal amplitude 
and the transverse amplitudes are suppressed by the factor mi^o/mg, which significantly 
conflicts with the LHCb and CDF observation f L (B° s K*°K*° , 00) ~ f ± (B° s K*°K*° , 00). 
Therefore, it is worthy to perform a detailed evaluation within QCDF, and check if the tension 
could be moderated by annihilation corrections. 

4.1 Within Scheme I 

With the annihilation parameters given by Eq. (1171) . the prediction of scenarios S4 for B® — > 
pp , K*°K*° and 00 decays are listed in the third column of Table HI which agree with the 



12 



Table 2: The numerical results for the branching fractions £>[xl0~ 6 ], the direct CP viola- 
tions AcpfxlO -2 ], longitudinal and transverse polarization fractions fi,± [xlO -2 ] for B® — > pp, 
K*°K*° and 00 decays in each scenarios. 
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former results of the QCDF [T5l [T7] . One may find that there are no significant direct CP 
asymmetries for these decay modes, and the longitudinal polarization fractions of B® — > pp 
decays are close to unity. The branching fraction of the B° s — > 00 decay agrees well with the 
experimental data. While, the default result B(B° S -> K*°K*°) ~ 11.0 x 10^ 6 is significantly 
smaller than LHCb measurement ~ 28.1 x 10~ 6 . However, one may notice that the uncertain- 
ties in the experimental measurement are still very large. For their polarization fractions, as 
expected above, the prediction of scenarios 5*4 implies —> K*°K*° ,00) ~ 0.71 ,0.56 > 

f±(B° ->■ K*°K*° ,00) ~ 0.13,0.21, which conflict with the LHCb and CDF observation 
f L (B° s ->■ K*°K*° , 00) « f L (B° s ->■ K*°K*° ,00). In the following, we would perform a de- 
tailed evaluations to check whether such a discrepancy could be moderated by the annihilation 
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Figure 4: The dependence of B{B° S -> K*°K*° , 00) and f L (B° -> K*°K*° , <jxj)) on the phases 
with different pY v values. The dashed lines correspond to the error bars (1.68a). 



corrections. 

The dependence of B(B° S -»■ K*°K*° ,00) and f L (B° s K*°K*° ,00) on the annihilation 
parameters is shown by Fig. HJ Comparing Fig. 0] (b) with (d), we find the constraint from 
fi,{Bs ~~ K*°K*°) on the annihilation parameters is weak due to its large experimental uncer- 
tainties. From Fig. [J] (d), one may find the phase 0^ y ~ —40° or 50° with any value of p v s v 
could be helpful to moderate the tension between the theoretical prediction and experimental 
result for }'l{.B s — > 00). While, with such cf^ v value, as shown in Figs@] (a) and (c), a small 
pY v ~ 1 is required by both B(B° S — > K*°K*°) and B(B° S — > 00). Furthermore, compared with 
such solutions, we also find a larger p v s v ~ 2 with a larger phase (f^ v ~ —150° or 160° are also 
possible solutions. A detailed numerical examination is performed in the following in due. 

Under the constraints from B(B° S -> K*°K*° , 00) and f L (B° s -> , 00), the allowed 



11 




(a) 



150 
100 
50 


-50 
100 
150 
0.0 





— ■ 5 VV D ; 


* 5 VV C 




<M' ^VVg 


5 VV A 



0.5 



1.0 



1.5 



2.0 



2.5 



(b) 



3.0 



Figure 5: The allowed regions for the annihilation parameters yy and pY v under the con- 
straints from B{B° S -)• K*°K*° , (jxf)), f L {B° s -)• K*°K*° , 00) (Fig. (a)) and their combination 
(Fig. (b)), respectively. 



regions for the annihilation parameters are shown in Fig. [5j As shown in Fig. [5] (a), the space 
of the annihilation parameters are strictly restricted. The upper and the lower limits of p^ v 
are dominated by B(B® ->■ 00) and B(B° -> K*°K*°), respectively. While, the ranges of vy 
are dominated by fiiB® — > 00). Finally, under their combined constraints, we could find four 
solutions shown by Fig. [5] (b). The corresponding numerical results are 



yv 



yv 



yv 



vv 



2.0 ±0.1 

1.0 ±0.1 

1.1 ±0.1 
2.0 ±0.1 



i^ = -154° ±4° (S VV A), 
^ = -36° ±9° (S VV B), 
i yy = 50°±7° (S VV C), 
i yy = 164°±4° (S VV D). 



(27) 



our predictions 
could reduced 



With the default values of the annihilation parameters given by Eq. (j2 
of scheme I are summarized in the Table [2j We find that fi{Bs ~^ K*°K*° 
to the experimental data by the annihilation contributions. Meanwhile, f±{B° s — y K*°K*° 
are also significantly enhanced, and agree well with the experimental data. However, similar to 
the case in scenarios S*4, the result of B(B® — y 00) is larger than the one of B(B° S — y K*°K*°) 
by a factor about 2. In the previous works, the theoretical predictions within both QCDF and 
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pQCD frameworks, for example 



^(4°^#) = (21.8l™)xl0- 6 , 
B(B° S -> K*°K*°) = (9.1±g;|±^8 3 ) x 10" 6 ; 

S(4°->00) = (16.7^ 1 8 3 )xlO- 6 , 

-)• = (6.6j££$) x 10- 6 ; 

B(B° -+ 00) = (35.3^^.7+0.0) x 10 _ 6 } 
B(B° S -»• K* i?* ) = (7.8li;^ji°;°) x 10' 



QCDF PS| (28) 
QCDF [IT] (29) 
pQCD [5J (30) 



present a similar result that B{B° S — > 00) > B{B Q S — >■ i^* ^* ), which is obviously inconsistent 
with the LHCb and CDF measurements Eqs. (J2SJ and ([26]) that £(5° -> 00) ~ B(B° S 
i^*°^*°). Such a theoretical situation could be easily understood from their amplitudes given 
by Eqs. flE} and (09]) in appendix A. The amplitudes of both B° s -> K*°K*° and 5° -»■ 00 
decays are dominated by the effective coefficients a^, and annihilation contributions to them 
are similar. However, an additional overall factor 2 is included in the amplitude of B® — > 00. 
So, 00) would be always larger than B(B° S -»■ K*°K*°). 

For the 5° — > pp decays, their longitudinal polarization fractions are always close to unity 
within the four solutions. While, the predictions of S VV B and S C for B(B° S — > pp) are 
significantly larger than the ones of S vv A and S VV D. So, the four solutions for the annihilation 
parameters given by Eq. fTSTl) could be tested by the upcoming LHC-b measurements of B(B® — > 
PP)- 

4.2 Within Scheme II 

With the infrared finite gluon propagator to deal with the endpoint divergences, the hard 
spectator and the annihilation corrections for B — > PP and PV decays have been evaluated in 
Ref. [Hj. While, the ones for B — > VV decays have not calculated until now. So, firstly, we 
would re-calculate these corrections within scheme II. With the same convention and notation 
as Refs. [Hj and [TBJ, the hard spectator scattering contributions can be expressed as 

Hr = VBf B tv f dxdyd&M 2 ) ^^^ , (31) 

m B m b F^ (0) Jo 1 ' + ^(q 2 ))xy 1 ; 
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for i = 1,2,3,4,9,10; 



rr- 2/b/v! t a A AC I 2^Bl<Pii.x)<p a2 {y) . . 



for i = 5,7; 



^ = f dxdvdZa a tf) t mM $tiM (33) 



for 2 = 6,8. In which, oo 2 (q 2 ) = M 2 (q 2 )/M B , q 2 = —Q 2 and Q 2 ~ —£xM B is the space-like gluon 
momentum square in the scattering kernels. The function $bi(£) is the B meson light-cone 
distribution amplitude, where £ is the light-cone momentum fraction of the light anti-quark in 
the B meson. In our following numerical evaluation, = N B C,(1 — £)exp — ^ M^-m b ) (£~ 

O?) 2 [23] is used. 

The longitudinal part of the annihilation amplitudes are given by 

/i j 
ofeofr/a s (g 2 ){$yi(£)$y2(2/) [7-= 2 , 2 x , ■ V1 — r + -7— 21 2\ i • d 
o I L (xy -o; 2 (g 2 ) +ie)(l - xy) y(:n/ - w 2 (g 2 ) + ze) J 

^^^ xy-^ + fe }' (34) 

^ 2 ° = « J tedy« a tf){* v M*vM [ (g _ _ u2{q2) x + m _ W) + x{x . _ J {q2) - . g) ] 

-r^^x)^) — 2 ^ j , (35) 

= tt/ dxdya s (q 2 )\r^<f> vl (x)$ V 2(y) 



(xy — uj 2 (q 2 ) + ie)(l — xy) 



^y^(y) {x ,_^ v +is){l _, y) }, (36) 



A f 3 '° = Ti f dxdya s (q 2 )^<^ vl (x)<l>v2(y)-7— 



2(1 + y) 



y(xy — uj 2 (q 2 ) + ie) 



x(xy — ui 2 (q 2 ) + ie) J ' 



(37) 
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and A{'° 2 = 0. The non-vanishing transverse amplitudes are 
A\~ = n 2mi ™ 2 I dxdya s (q 2 )l(j) bl (x)(j) b 2(y)[- 



m 2 B J o I L (l — xy)(xy — u 2 (q 2 ) + ie) 



2 2 

xy xy 



(1 — xy)(xy — u 2 (q 2 ) + is) 2 (1 — xy) 2 (xy — u 2 (q 2 ) + ie) 
' (xy-u 2 (q 2 )+ie) 2 ^ } ' (38) 



A\ = n 2mi ™ 2 I dxdya s (q 2 )((j) al (x)(f) a 2(y)[- 



m 2 B J I L (l - xy)(xy - u 2 (q 2 ) + ie) 



2 2 

x y x y 



(1 — xy)(xy — u 2 (q 2 ) + ie) 2 (1 — xy) 2 (xy — ui 2 (q 2 ) + ie) 

r] } , (39) 



(xy — u 2 (q 2 ) + ie)' 



A l i = tt I dxdya s (q 2 )\^-r^ 2 (f) al (x)(t>2(y)- 



o 



m 2 x (xy — uj 2 (q 2 ) + ie)(l — xy) 



r?<ti(x)Mv) ,. 2 ,„ 2 x I 737i ^ > (40) 



mi x (xy — cu 2 (q 2 ) + ie)(l — xy) 



A3 = Ti dxdya s (q 2 ) 
J 



2mi v 2 > / \ / 1 / \ 



-r 2 a i(x)0 2 (y)- 



m 2 x y(xy — u 2 (q 2 ) + ie) 



+ ^l r v^± {x)My) 1 — — } , (41) 

mi x(xy — <jj l (q l ) + is) > 

= 7T g / GbA/a S (g 2 W 4> a l(x)4>a2(y) [ 



+ 



J o I L (l — xy)(xy — u 2 (q 2 ) + ie) 

xyy xyy 



(1 — xy)(xy — oo 2 (q 2 ) + ie) 2 (1 — xy) 2 (xy — u 2 (q 2 ) + ie) 
1 



y 2 (xy — co 2 (q 2 ) + ie) y(xy — ui 2 (q 2 ) + ie) 



] } , (42) 



= 71 2 — / dxdya s (q 2 )<(j) bl (x)(j) b 2(y)[ 



+ 



m% J o I L (l - xy)(xy - w 2 (g 2 ) + ie) 

xxy xxy 



(1 — xy)(xy — u 2 (q 2 ) + ze) 2 (1 — xy) 2 (xy — u 2 (q 2 ) + ie) 
1 xy 



] } , (43) 



x 2 (xy — ui 2 (q 2 ) + ie) x(xy — oo 2 (q 2 ) + ie) 2 - 
where q 2 ~ xyM B is the time-like gluon momentum square. As found in Ref. |14j . the hard- 
spectator scattering contributions are real, but the annihilation contributions are complex with 
a large imaginary part. Their contributions are dominated by the value of the effective dynam- 
ical gluon mass scale m g . 

With the default values of the input parameters, the dependence of B(B° S — > K*°K*° , 0</>) 
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Figure 6: The dependence of B(B° S K*°K*° ,00) and f L ,±(B° -> K*°K*° ,cjxf)) on the 
effective dynamical gluon mass scale. The dashed lines correspond to the error bars (1.68<t). 

and /x — > j^* ^* ; 00) on the parameter m g is shown by Fig. From Figs. M (b) and (c), 
we find that the QCDF predictions for fn±.)(Bs —> K*°K*° ,00) could be reduced (enhanced) 
to the experimental measurements with a small effective dynamical gluon mass m g < 0.4GeV, 
which leads to a large annihilation contribution. Meanwhile, as Fig. |6] (a) shows, a small m g 
is also allowed by the constraint from B{B° S — > K*°K*°). However, such small m g value would 
induce too large B(B° S — > <fr<f)), which is much larger than the experimental data. 

With a default m g value 0.5GeV, we present our predictions for the observables in the last 
column of Table [3 We find that our prediction B(B° S -»• K*°K*°) = (20.6±f;|) x 10~ 6 agrees 
well with the experimental data (28.1 ± 6.5) x 10 -6 . However, unfortunately, B(B® — > <fxj)) is 
enhanced much by the annihilation corrections, which lower limit (49.9 — 13.3) x 10~ 6 is about 
1.6a larger than the CDF measurement (23.2 ± 8.4) x 10 -6 . In fact, with any m g value, as 
analysis in the last section, the theoretical prediction of B(B° S 00) is always much larger 
than B(B° S — > A'* ^* ), which is also can be seen from Fig. [6] (a). Because a small m g , which 
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is help to accommodate the discrepency of fi,±(B® — > K*°K*° , 00) between the theoretical 
predictions and experimental data as Figs. [6] (b) and (c) showing, is excluded by the constraint 
from B(B° S — > 00), the predictions of scheme II for fi(±){Bs — > K*°K*° , 00) are larger (smaller) 
than the experimental measurements, which is similar to the results of scenarios 54. For the 
B® — > pp decays, their branching fractions in scheme II are significantly larger than the ones in 
scheme I. So, the upcoming measurements on B(B® — > pp) could give a judgment on the two 
schemes. 

5 Conclusion 

Motivated by the most recently observed large branching fraction of the pure annihilation decay 
B° s — > 7r + 7r~ and large transverse polarization fractions in the B° s — > K*°K*° , 00 decays, we 
revisit the hard spectator and annihilation corrections in non-leptonic B° s decays within QCD 
factorization approach. In this paper, two schemes (parameterization and using an infrared 
finite gluon propagator) to model the effects of the end-point singularity in hard spectator and 
annihilation corrections are evaluated. In our numerical evaluations, all of the pure annihilation 
and well measured B® decays are studied in detail simultaneously. Our main conclusions are 
summarized as: 

• For B° s — y PP decays, within scheme I, due to the large B(B° S — > ti + tx~) measured by 
CDF and LHCb collaborations, a large annihilation parameter p p p is required. Under the 
constraints from B(B° S — > 7r + 7r~ , tt~K + , K~K + ) and Acp(B° s —> n~K + ), two solutions 
pP p 2.5(3.5) and 0f p ~ —84° (116°) are found, which are significantly larger than 
these values p pp ~ 1 and PP ~ —55° adapted in the literature. With the obtained two 
solutions of the annihilation parameters given in Eq. (Tl8|) . all of the the QCDF predictions 
for B° s — > PP decays agree well with the available experimental data. 

• For B® — > VV decays, within scheme I, the measured small fi{Bg — > K*°K*° ,00) and 
large f±{B° s — > K*°K*° ,00) could be accommodated by the annihilation contributions. 
With the constraints from B(B° S K*°K*° , 00) and /z,±0B° -> K*°K*° , 00), we find 
four solutions of the annihilation parameters p^ v and (f^ v given by Eq. (|2T|) . Some of 
these solutions will be confirmed or refuted by the upcoming LHCb measurement on 
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B® — > pp decays. 

• Within scheme II, using the effective gluon mass m g = 0.48GeV, QCDF predictions for 

— > PP decays are found to be in good agreement with the available experimental 
results. Furthermore, some of its predictions are different from the ones in scheme I, such 
as the branching fractions and direct CP asymmetries of B° s — > pir decays, which will be 
judged by the upcoming LHCb and proposed super-B experiments. For the B° s — > VV de- 
cays, Jl,±{Bs ~~ >* K*°K*° , 00) could be accommodated by the annihilation contributions 
with a small m g < 0.4GeV, which unfortunately is excluded by B(B° S — > 00). 

• Within both scheme I and II, B(B° S — > 00) is always larger than B(B® — > K*°K*°), which 
significantly conflicts with the LHCb and CDF observation B(B° S — > 00) ~ B(B° S — > 
K*°K*°). A similar situation is also presented in pQCD approach as summarized in 
Ref.pl] and Eqs. (TJHl ES [30]). Thus, the present experimental result B(B° 00) ~ 
B{B° S — > x*°K*°) rises a challenge to theoretical approaches for B nonleptonic decay. 
The further refined measurements and theoretical studies are required to resolve such a 
possible anomaly. 
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AppendixA: Decay amplitudes with QCDF 



For selfconsitence of this paper, the following decay amplitudes are recapitulated from Refs. [12 
[15]. The amplitudes of B s — > up decays are 



B-u-p 
B 



°4 2 4,EW 



B 



pTT 



Spuh + l?A + b p A 



4,EW 



ix p 



5 pu bi + 2If 4 + -^ iEW 



B 



1 



Spuh + 2bl + -bl 



4,EW 



(44) 
(45) 
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The amplitude of B s — > n~ p + decay is obtained from the first expression by interchanging 
(it) -h- (p) everywhere. The expressions for the B s — > Tin and pp amplitudes are obtained by 
setting (p) — > (7r) and (n) — > (p), respectively 

The decay amplitudes of B s — > n~K + , K~K + , and K*°K*° , (fxp decays are : 



A 



B s ^n~K+ 



A 



B a ^K~K+ 



A 



Akit 

Bk-k+ 
+A K+K 

-\-A k *ok*o 



<5piA + &4 + ^4 EW 



2^3,EW 



4,EW 



2 ^.EW 

1 



ryP - -ry p + R p + B p --B p - -B p 



.4 



2A 



«3 + (X A 



2"3,EW 2" 4 ' EW ^3 + ^4 2^3,EW 2^,EW 



-O 



-a 



p 

4,EW 



(46) 
(47) 

(48) 
(49) 



The explicit expressions for the coefficients = af (iWiM 2 ) and /3f = (3f (M X M 2 ) can be found 
in Ref. 



AppendixB: Theoretical input parameters 

For the CKM matrix elements, we adopt the CKMfitter Group's fitting results [25] 



p = 0.144 ± 0.025 ,rj = 0.342lg;g^ , A = 0.812+^ , A = 0.22543 ± 0.00077 . 



0.013 



(50) 



As for the quark mass, we take [2nj 123 [2S] 



m u = ma = m s = 0, m c = 1.611q 12 GeV, 



m b = 4.78l{j$ GeV, m t = 172.4 ± 1.22 GeV. 



(51) 



for the pole masses and 
m s (p) 



m q (n) 

m b (m b ) 



27 A ± 0.4 , m s {2 GeV) = 87 ± 6 MeV, m c (m c ) 
4.19i^ GeV , m t (m t ) = 164.8 ± 1.2 GeV , 



1.27^ GeV 



(52) 



for the running masses, where m„ = m Ul 
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The decay constants are [221 [2§J [3U] 

f Bs = (231 ±15) MeVJ Bd = (190 ± 13) MeV , 

f n = (130.4 ± 0.2) MeV , f K = (156.1 ± 0.8) MeV , f p = (209 ± 2) MeV . (53) 

As for the B-meson lifetimes, tb a = 1.525 ps [26] and tb b = 1.472 ps [26] are used. Furthermore, 
the form factor F B ^(0) = 0.258 ±0.031 [3U] is used to evaluate the amplitude of B° d 7i + K~ 
decay. For B Q S n~K + , K~ K + decays, we shall use Fq s ^ k (0) = 0.24 obtained by both lattice 
and pQCD calculations and suggested by Ref. [17]. For the other decay constants and form 
factors related to B° s — > K*°K*° and B° s — > 00, we choose the similar inputs used in Ref. |15j . 
In which, these values follow the QCD sum rule calculation Ref. [30], but some modifications 
within theoretical errors are made to improve the description of data. Their values are 

f K * = (218 ± 4) MeV , /^»(2GeV) = (175 ± 25) MeV , 

U = (221 ± 3) MeV , /^(2GeV) = (175 ± 25) MeV , 

a b s ^k* = Q 38 ± q Q3 ; F Bs^i<* = 0.53 ± o.05 , F^ R ' = 0.00 ± 0.06 

= 0.38±g;^ , F* ^ = 0.65±g;^ , F^ = 0.00 ± 0.06 . (54) 
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